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Q> ; Abstract 

^ ■ A theory of electron counting statistics in quantum transport is presented. It 

Ch , involves an idealized scheme of current measurement using a spin 1/2 coupled 

O 

to the current so that it precesses at the rate proportional to the current. 
Within such an approach, counting charge without breaking the circuit is 
possible. As an application, we derive the counting statistics in a single chan- 
nel conductor at finite temperature and bias. For a perfectly transmitting 
channel the counting distribution is gaussian, both for zero-point fluctuations 
and at finite temperature. At constant bias and low temperature the distribu- 
tion is binomial, i.e., it arises from Bernoulli statistics. Another application 
considered is the noise due to short current pulses that involve few electrons. 
We find the time-dependence of the driving potential that produces coherent 
noise-minimizing current pulses, and display analogies of such current states 



with quantum-mechanical coherent states. 
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I. INTRODUCTION 

Quantum transport in nanostructures has been a subject of many recent studies 0. 
Transport properties like Ohmic conductivity can be understood in terms of the quantum 
scattering problem in the conductor, which provides a theory of quantum coherence of 



transport [Q. Fluctuations of electric current due to the discreteness of electric charge 
are intrinsic to quantum transport It has been found that current fluctuations have 

interesting properties reflecting profound aspects of underlying quantum dynamics l^j-Q. 
For example, the quantum noise caused by a dc current is reduced below classical shot noise 
level 0-0. This suppression has been understood as an effect of enhanced regularity of 
transmission events due to Fermi statistics [0. Besides theoretical interest, such phenomena 
may lead to applications. Given the development of nano-technologies, the transmission of 
signals by single- or few-electron pulses will become common, and then one will see the 
quantum statistics of current working. 



In this paper we update the theory of quantum measurement of electric current |]TT|. Our 
goal is a complete description of charge fluctuations, rather than developing measurement 
theory (see Sees. II, III). We shall derive a microscopic formula for electron counting distri- 
bution (see Sec. Ill, Eq. (0), and Sec. IV, Eq. (|26|) ) that can be used for any system, e.g., 
with an interaction or with a time-dependent potential [jl2|. As an application, we test the 



method on the statistics in a single channel ideal conductor for non-equilibrium and equilib- 
rium noise at finite temperature, and for zero-point equilibrium fluctuations (Sees. IV, V). 
In particular, the fluctuations of a dc current at zero temperature are found to be binomial 
(Sec. VI) with the probabilities of outcomes related with transmission coefficients of elastic 
scattering in the system, and with the number of attempts = eVt/h, where V is applied 
voltage, and t is the time of measurement. 

Another property of quantum noise that does not have classical analog is its phase 
sensitivity OJT^. For the current correlator ((j(ti)j(t2)))+ it results in a periodic sinusoidal 



dependence on Faraday's flux due to applied voltage, $ = c fl^ V{t)dt, with the period 



$0 = hc/e. The phase sensitivity manifests in singularities of the low frequency noise power 
in a junction driven by ac and dc signals together [|T5|. 

Even more remarkable is the behavior of current fluctuations due to a pulse of voltage 
T3| , p!6|| . Total charge that flows through the conductor due to a voltage pulse fluctuates in 



such a way that the mean square fluctuation diverges whenever the flux of the pulse is not 
an integer: = | V{t)dt ^ 2?™. On the other hand, for = 27m the fluctuation of 
the transmitted charge is finite (Sec. VII). This result has simple interpretation in terms of 
the Anderson orthogonality catastrophe theory, since the flux (p enters the time dependent 
scattering matrix of the conductor through the forward scattering amplitude. 

With this, one is led to address the issue of current states that minimize the current 
fluctuations at fixed mean transmitted charge 117, TB|. It is found in Sec. VIII that such 



states are produced by time-dependent voltage of the form 

a sum of Lorentzian pulses of unit flux each. It is remarkable that the minimal noise 
due to such sequence of pulses is independent of the pulse positions tk and widths Tk, 
which leads to obvious parallels with solitons in the theory of non-linear integrable systems. 
The noise minimizing current states can be compared to the coherent states that minimize 
the quantum-mechanical uncertainty. Apart from obvious similarity, there is a difference: 
the coherent current states are many-body time-dependent scattering states. Their role in 
transport is an interesting subject of future work: one expects that representing many-body 
states as a superposition of these coherent states has an advantage similar to that provided 
by coherent states of one particle. 



II. MEASURING ELECTRIC CURRENT 



Instantaneous measurement is described in quantum mechanics by wavepacket reduction 
that involves projecting on eigenstates of an observable, represented by a hermitian operator. 
A different kind of measurement, extended in the time domain, is realized in detectors and 
other counting devices. It is known that in such cases a certain revision of the measurement 



description is necessary. A famous example is the theory of photon detectors |T8[ in quantum 
optics. Due to Bose statistics, photons entering a photo-counter are correlated in time, and 
this makes the theory of photon detection a problem of many-particle statistics. For a single 
normal mode of radiation field the probability Pm to count m photons over time t is given 
by 

Prn = ^^(:(a+a)™e-''*'^^'^:), (2) 

where and a are Bose operators of the mode, rj is the counter efficiency parameter, and 
(...) stands for the average over a quantum state. The normal ordering : . . . : is an important 
element of the formalism. Physically, it means that, after having been detected, each photon 
is destroyed, e.g., it is absorbed in the detector. Instead of the probabilities (0), it is more 
convenient to deal with the generating function 

X(A) = Pmc''"^ ■ (3) 

m 

For the single normal mode Eq.(^ leads to 

X(A) = ( : exprit{e'^ - l)a+a : ) . (4) 

Eqs.(||,|,H) account very well for numerous experimental situations ||T9|. Particularly inter- 
esting is the case of a coherent state |z), a\z) = z |z), where z is a complex number. It 



corresponds to the radiation field of an ideal laser, and with Eq.(^) one easily gets Poisson 
counting distribution, 

Pn. = ^-^e-'\ J = V\^\\ (5) 
ml 

which describes the so-called minimally bunched light source. 

Similar to the photon detection, electric measurement is performed on a system contain- 
ing an enormous number of particles — in this case fermions — and thus one expects the 
effects of Fermi statistics to be important. Also, the duration of electric measurement is 
typically much longer than the time it takes the system to transmit one electron by micro- 
scopic tunneling, scattering, or diffusion. Apart from these similarities, there is, however, 
a crucial difference from the photon counting: the number of electrons is not changed by 
the current measurement, since electric charge is conserved. This has to be contrasted with 
absorption of photons in photo-counters. Related to this, there is another important dif- 
ference: at every detection of a photon, its energy hu is taken from the radiation field, 
which makes plain photodetectors insensitive to zero-point fluctuations of electromagnetic 
field. On the contrary, the measurement of current fluctuation is usually performed without 
changing energy of the system, which makes the zero-point noise an unavoidable component 



of any electric measurement ||20|. (Let us emphasize that the difference has nothing to do 



with the type of quantum statistics, Fermi or Bose. Rather it is the difference between the 



two kinds of measurement, e.g., see where counting of fermions was discussed using an 
optical-like counter that has to capture an electron in order to detect it.) 

In the classical picture, the measurement gives the charge Q{t) = jQj{t')dt' transmitted 
during the measurement time t. The probabilities Pm of counting m electrons can then be 



obtained by averaging S{Q{t) — me) over the state of the system. In a quantum problem 
electric current is an operator, and since currents at different moments do not commute, 
the operator of transmitted charge Q{t) = jQj{t')dt' generally does not make any sense. 
Instead, since we are interested in higher-order statistics of current fluctuations, beyond 
(j(t)) and {{j(ti)j(t2)))+, in order to compute electron counting ditribution, we have to 
include the measuring system in the quantum Hamiltonian. Our approach is motivated 
by the example of the quantum mechanical systems with strong coupling to macroscopic 
environment, introduced by Leggett, that can be treated consistently only by adding the 
"measuring environment" to the quantum problem . 



For that we introduce a model quantum galvanometer, a spin 1/2 that precesses in the 
magnetic field B of the current. For a classical system, the rate of precession is proportional 
to B{t), and B{t) is proportional to the current I{t): B{t) = const I{t). Therefore, the 
precession angle of the spin directly measures transmitted charge 5Q = jQl{t')dt'. We adopt 
the same measurement procedure for the quantum circuit, i.e., we include in the electron 
Hamiltonian the vector potential due to the spin: 

A(r) = -fiB X , (6) 
\r\ 

where a = [a^, Cy, a^) are Pauli matrices. Thus we obtain a Hamiltonian describing motion 
of electrons, the measuring spin, and their coupling. Now, according to what has been 
said, we have to solve dynamics of the spin in the presence of the fluctuating current, find 
the distribution of precession angles, and then interpret it as a distribution of transmitted 
charge. Of course, a question remains about the back effect of the spin on the system, as in 



any other problem of quantum measurement. However, as we find below in (^) and ( |19| ) 



only the phase of an electron state is affected by the presence of the spin, not the amplitude. 
Moreover, the phase will change only for the transmitted, but not for the reflected wave. 
As a result, the probabilities we obtain do not depend on the coupling constant of the spin. 
This justifies the assumption that the spin measures charge transfer in a non-invasive way. 



It is worth remarking that our scheme resembles the "Larmor clock" approach ||23[ to 
the problem of traversal time for motion through a classically forbidden region. In this 
problem one is interested, e.g., in the time spent by a particle tunneling through a barrier. 
The Larmor clock approach involves an auxiliary constant magnetic field B added in the 
classically forbidden region, and a spin 1/2 carried by the particle that interacts with the 
field: Tiint = —<^zB. The precesson angle of the spin measures traversal time. Comparing 
the two approaches is very tutorial: see Appendix A, where the Larmor clock is reviewed. 



III. SPIN 1/2 AS A GALVANOMETER 

Having clarified our motivation, we proceed semi-phenomenologically and choose a new 
vector potential in the spin-current interaction —\jA. We replace the Ampere's long-range 
form d^) by a model vector potential 

A,{r) = ^d^V,e{f{r)-h) (7) 

concentrated on some surface S defined by the equation /(r) = /q. Here $0 = hcje, A is a 
coupling constant, /(r) is an arbitrary function, and, as usual, the step-function 6{x) = 1 for 
X > 0, for x < 0. The surface 5* defines a section of the conductor on which the interaction 
is localized: 

/I Xh 
— jA d?r = -—aJs , (8) 
c 2e 

Q 



where Is = Jsjds, i.e., the spin now is coupled to the total current through the section 
S. With the choice (|^ of the vector potential one can study current fluctuations in an 
arbitrary section of the conductor. Another advantage of the phenomenological Eq.(|^ is 
that it involves only one Pauli matrix, which makes the spin dynamics essentially trivial. 
The choice of the quantization axis of the spin is arbitrary since (|^) will be the only spin- 
dependent part of the Hamiltonian. Finally, another advantage of the form (H) is that by 
switching from the smooth function (|]) to the singular form (|^) we enforce integer values of 
counted charge. To understand this, let us note that in the "fuzzy" case (|^) the measurement 
can start at the moment when one of the electrons is located somewhere in the middle of 
the volume where A ^ 0, and then a fractional part of electron charge will be counted. On 
the contrary, in the "sharp" case (|^, the spin responds to the presence of an electron only 
when it crosses the section S. We shall see below in a microscopic calculation that integer 
values of charge follow automatically from gauge invariance, since the form (|^) is a gradient 
of a scalar. 

Thus we come to the Hamiltonian 

= H{p, r), p^=p^- -Ai , (9) 

c 

— * 

where the spin- dependent A is taken in the form (|^). An essential feature of our approach is 
that we treat the constant A of coupling between the spin and the current as a variable, i.e., 
we consider the spin precession as a function of the parameter A. The reason is that, unlike 
the photon counting problem, our measurement scheme directly generates the function x(A), 
and then the counting probabilities Pm are obtained by reading Eq.(Q) backwards. 

At this point we are able to formulate our main result. Let us define a new Hamiltonian 

in 



Hx = niP,r), P.=Pi- \\h V.eifir) - /o) , (10) 

simply by supressing in Eq.(|^). The Hamiltonian H\ involves only quantities of the 
electron subsystem. Below we show that by measuring precession of the spin coupled to the 
current, one obtains the quantity 

;^(A) = (e^W-A*e-^«A* ) . (11) 

Here the brackets (...) stand for averaging over initial state of electrons. Note that xW is 
written in terms of a purely electron problem, not involving spin variables. We shall find 
that the function xW defines the result of any measurement of the spin polarization after 
the time t when the spin-current coupling is turned off. Moreover, we shall see that the 
function (|Tl|) has the meaning of a generating function of electron counting distribution. 



i.e., the Fourier transform of xW gives counting probabilities, entirely analogous to (H). 

Our goal now will be to express evolution of the spin in terms of quantities corresponding 
to the electron system. The interaction is given by Eqs.(0),(|^). Suppose that the measure- 
ment started at the moment and stopped at the time t, i.e., the spin-current interaction 
is turned on during the time interval < r < t. Let us evaluate the density matrix Ps{t) of 
the spin, right after it is disconnected from the circuit. We have 

Ps{t) = tre(e-*^-*pe^«-*) , (12) 

where p is initial density matrix Pe®Ps at t = 0, Pe is initial density matrix of electrons, and 
tre(...) means partial trace taken over electron states, the spin indices left free. In terms of 
the spin variables, the operator e"*^''* is a function only of a^, and hence it is diagonal in 
spin: (T|e~*''^''*|i) = (i|e^*''^'^*|t) = 0. In other words, if initially the spin is in a pure state, 
up or down, it will not precess. For Ps{t) this remark yields: 

1 1 



Ps{t) 



(13) 



ptt(o) x(A)Pn(o)' 
xi-m^io) Pu(0) . 

Here xW = tre(e~*^^*pee*^"^*), where e~*^^* is the evolution operator for the problem (|T0|). 
Now, after the spin degrees of freedom are taken care of by (p!3D , we are left with a purely 
electron problem, that involves only electron degrees of freedom, but not the spin. By using 
the cyclic property of the trace tre (...), one can show that xW iii Eq.(|T3D is identical to 



11 



In principle, any entry of a density matrix can be measured, and therefore the quantity 
x(A) is also measurable. In order to make clear the relation of xW with the distribution of 
precession angles, let us recall the transformation rule for the spin 1/2 density matrix under 
rotation by an angle 6 around the axis: 



Ptt Pn 



nid: 



e Pit PU 
By combining this with Eq.(|D we write Ps(t) as 

Ps(t) = ^PmTle=m\{p) 



(14) 



(15) 



which assigns to Pm the meaning of the probability to observe precession angle m\. Let us 
finally note that such interpretation of is consistent with what one expects on classical 
grounds, because for a classical magnetic moment a interacting with the current according 
to (0), the angle 9 = \ corresponds to the precession due to a current pulse carrying the 
charge of one electron. 
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IV. SINGLE-CHANNEL CONDUCTOR. GENERAL FORMALISM 

In order to see Eq.(|lT]) working, let us consider an ideal single channel conductor, i.e., 
the Schrodinger equation 



dip 



(16) 



in one dimension, where the potential U{x) represents scattering region and the vector 
potential is inserted according to (|^ at the x = section. In order to describe transport, we 
shall use scattering states, left and right. Their energy distributions nL{R){E) are equilibrium 
Fermi functions with temperature T and chemical potentials shifted by eV , — = eV, 
representing a dc voltage. 

For the problem one can write time dependent scattering states as 

gifcx _^ 5^e~*''^, X < —a/2 
X > a/2 



-iEkt 



ikx 



i'RA^.t) = e 



^-i\/2^^^-ikx^ X < -a/2 



(17) 



g-ifcz _^ iJ^e'''^, X > a/2 
where a is the width of the barrier, and Al^r and Blr are the transmission and reflection 
amplitudes in the absence of the spin vector potential. To make expressions less heavy, we 
supress electron spin. The phase factors e^*^/^ in ( [T7| ) are found immediately by observing 
that the vector potential in the Schrodinger equation can be eliminated by the gauge trans- 
formation iIj{x) — >■ exp(iA/2 6{x))iIj{x). The scattering amplitudes form a unitary matrix: 

We will study the range of small T, eV <C Ep, when only the states close to the Fermi 



(18) 



level are important. In this case, there is an additional simplification because the states 



15? 



near Fermi energy have almost linear dispersion, and thus all wavepackets travel with the 
speed vp without changing shape. Then, following Landauer and Martin [p, instead of the 
usual scattering states (|T^, it is convenient to use their Fourier transform. By ignoring the 
energy-dependence of A^^r and -Bl,_r, which is equivalent to saying that the scattering time is 
negligible, and assuming that the dispersion is stricktly linear, one obtains the representation 
of scattering in terms of time- dependent scattering wave packets 



5{X-), t < T 



t < T 

V'i?,r(x,t) = <! (19) 

e"^/^AR5{x+) + Br5{x^), t > t 
where x± = x ± Vpit — t). Here r is the packet arrival moment, at which it is scattered. It 



is straightforward to verify orthogonality of the states (|T^). Any time- dependent electron 



state can be written as a superposition of the states ([T9|), with the arrival time r serving in 



this representation as a label in the continuum of states, like k in (0). The assumption that 
the scattering amplitudes are energy-independent (and thus the scattering takes no time) is 
equivalent to replacing the barrier U (x) of finite width by Uo6{x) and is consistent with the 
closeness of relevant energies to Ep. 



Second-quantized, electron states (|TUp lead to ip{x,t) = ■ipL{x,t) + ijjji{x,t) with 

'>pL{R) (X, t) = J2'^L{R),r (x, t)ci{2),r , (20) 



where ci^r and C2,r are canonical Fermi operators corresponding to the states (|19|) , the left 
and the right respectively. One checks that fermionic commutation relations for Ci(2),r; 

4,rCj,T' + Cj, r'C+^ = 5ij5 (r - r') (21) 

14 



Ci,TCj,r' + CjyCi^r = , C+^C+^' + Str'^^r = ' (22) 

yield the usual commutation relations for il'L(R){x,t). From that one finds the meaning of 
the summation in pO|): = /^•••'^t- Mathematically, in this paragraph we defined 

T 

second-quantized ipi^) in (P^- 

The advantage of introducing the basis of the wavepackets (0),(^) is that now it is 
straightforward to write the many-particle evolution operator through the single-particle 
scattering matrix S\: 



exp y dTY^\\i[Sx]ijcj.^Cj^r , (23) 



where S\ is given by (]TB|). To verify (PBD, let us note that in the wavepacket representation 
(PD|), according to Eqs. (PT|), Fermi correlations occur only for the pairs of left and right states 
that scatter at the same instant of time. For each of such pairs the evolution operator e"*"^^* 
is 1 if both states are occupied or both are empty, otherwise it is given by the single-particle 
scattering matrix (0). 

Using similar arguments, we compute 



exp f drY, WijC+^Cj^r , (24) 



where = S_lS\ is readily obtained from (|I8D : 



(25) 



2t sin XBrAl e-'^\AR\' + \BR\ 
Using unitarity of and commutation rules for Ca,T one can rewrite ( ^4[) in terms of normal 
ordering: 



t 

^in^,t^-m,t = . exp / cir^[e^ - l]ijC+^c,-, : (26) 



^3 
IK 



This form is ready to be plugged into Eq.([Tl|) and averaged over the initial state. Let us note 
the striking similarity of two formulas obtained by different means: the fermionic Eq.(p6D 
and the bosonic Eq.(^. 

Also, let us mention that the periodicity of the matrix ( ^51) in A ensures periodicity of 
x(A), and thus guarantees integer values of charge. 

V. EQUILIBRIUM FLUCTUATIONS 

Let us start with a simple one-particle example. Consider a particle in the state c^^|vac), 



which corresponds to scattering at the moment r. In this case, from (pOf ) and (|Tl|) one gets 

'e*^U|2 + |5|2 forO<r<t: 



X(A) = 



(27) 

1 otherwise; 



\A\ = \Al\ = \Aji\, \B\ = \Bl\ = \Br\. Evidently, according to Eq.(|), this simply means 
that for the scattering occurring during operation of the detector, the counting probabilities 
are identical to the one-particle scattering probabilities, as it should be expected. 

Now, we consider current fluctuations in an equilibrium Fermi gas. First, let us assume 
perfect transmission: Bl{r) = 0. Then Eq. fpS]) gives W = iXaz, and thus Eq.(p^ becomes 



expiX f {cf^ci^r - ctrC2,T)dT , (28) 

JO 



i.e., the right and the left states separate. 

We observe that averaging of ( PBD over the Fermi ground state is identical to that per- 
formed in the orthogonality catastrophe calculation P^ , P3[] . Hence, averaging of (^) can be 
done by using the bosonization method [^] that replaces the fermionic Hamiltonian by a 
bosonic one. (The calculation is described in Appendix B.) In the bosonized representation 
one has to do a simple gaussian average, which gives 



X(A) = e-^^/(*'^) , (29) 
where A/27r + 1/2 = [A/27r + 1/2J, with [...J being the fractional part. The function 

I Tt//i- 2^1n27rT(5 , t^h/T 
where 5 is an ultraviolet cutoff time, of the order of h/Ep- At long times, according to @, 

this leads to gaussian counting statistics. 

As a side remark, the distribution given by Eq.(p9D also gives a solution to another 
problem: the statistics of the number of fermions inside a segment of fixed length in one 
dimension. The relation is immediately obvious after one assigns to r in Eq.(p8D the meaning 
of a coordinate on a line. Thus, in this problem the statistics are gaussian as well. 

Now, it turns out that the general case of non- vanishing reflection, B ^ 0, can be reduced 
to (^) by a canonical Bogoliubov transformation of c^^r making the quadratic form in (plf ) 
diagonal. The transformation is related in the usual way with the eigenvectors of the matrix 
W. Thus, we come to Eqs.(p8D,(p9D with A replaced by A,,: 

sin — = 1^41 sin — . (32) 
2 ' ' 2 ^ ^ 

The counting statistics in this case are non-gaussian: 

X(A) = e-^^/(*'^) . (33) 
One checks that the second moment of the distribution 



agrees with the Johnson-Nyquist formula for the equilibrium noise 

17 



= 2\A\'f{t,T) (34) 

A=0 



VI. STATISTICS OF A DC CURRENT: QUANTUM SHOT NOISE 

Let us consider non-equilibrium noise. In this case, due to the asymmetry in the distribu- 
tions, nL(R){E) — {exp{E±^eV) /T + generally one cannot uncouple the two channels 
by a canonical transformation. We calculate the statistics within an approximation that 
ignores the effect of switching at r = and r — t. Let us close the axis r into a circle of 
length t, which amounts to restricting on periodic states: 

^(r) =-0(T±i) . (35) 
For the periodic problem, by going to the Fourier space, one has 



X(A) = n [l + \A\\e-'^ - l)n4E,)(l - n«(Efc)) 
kez 

+\A\\e'^-l)nR{Ek)il-nLiEk))] , (36) 



where E^ = 27Thk/t, k is an integer. For large t, t ^ fi/T or t 3> h/eV, the product is 
converted to an integral: 

ln(x(A)) = dE In (l + \A\\e-'^ - 1) 

xni(l - nn) + {A^e'^ - l)n^(l - n^)) . (37) 

We evaluate it analytically, and get 

x(A) = exp (— tTM+M_//i) , (38) 

where 

u± — V ±cosh~^(|A|2cosh(i; + a) + l^l^cosh^;) , (39) 



IS 



V = eV/2T. The answer simplifies in the two hmits: T ^ eV and eV ^ T. In the first case 
we return to the equihbrium result (|33|). In the second case, corresponding to the recently 
discussed quantum shot noise we have 

X(A) = {e''^\A\^ + \B\Y^\'/\ e = sgnV , (40) 

Analyzed according to Eq.(|^), this xW leads to the binomial distribution 



p = |v4p, q = l-Bp, N = e\V\t/h. One checks that the moments (m) = pN and ((wi^)) = pqN 
correspond directly to the Landauer formula and to the formula for the intensity of the 
quantum shot noise 0-|5|. The correction to the statistics due to the switching effects is 
insignificant [|T^ . 



VII. NOISE DUE TO A VOLTAGE PULSE: ORTHOGONALITY CATASTROPHE 

Here we consider the fluctuations of current in a single- channel conductor induced by a 
voltage pulse. The result will be that the dependence of the fluctuations on Faraday's flux 
$ = — c / V{t)dt contains a logarithmically divergent term periodic in $ with the period 
$0 = hc/e. The fluctuation is smallest near $ = n^Q. The divergence is explained by 
a comparison with the orthogonality catastrophe problem. The $o "Periodicity is related 
with the discreteness of "attempts" in the binomial statistics picture of charge fluctuations 
presented above. 

Initially, the orthogonality catastrophe problem emerged from the observation that the 
ground state of a Fermi system with a localized perturbation is orthogonal to the non- 

1Q 



perturbed ground state, no matter how weak the perturbation |2^. Originally, the discus- 
sion was focused on the purely static effect of Fermi correlations on the ground state that 
leads to the orthogonality, but then it shifted to dynamical effects. When a sudden localized 
perturbation is turned on in a Fermi gas, the number of excited particle-hole pairs detected 
over a large time interval t diverges as Int/r, where r is the time of switching of the per- 
turbation. This effect leads to power law singularities in transition rates involving collective 



response of fermions, such as X-ray absorption in metals p8| , p4| . In this section we present 
an application of the orthogonality catastrophe picture to the electric current noise. 

Let us consider a single channel conductor in an external field described by the one- 
dimensional Schrodinger equation 

d ^ 

i—i){x,t) = nij{x,t) , 

T^=l[-^§^--Mx,t)J + Uix), (41) 

where the potential U{x) represents the scattering region and A{x,t) is the vector potential 
corresponding to the applied pulse of electric field. Since the pulse duration r is assumed 
to be much longer than the time of scattering, one can treat the vector potential as static 
and apply a gauge transformation in order to accumulate the flux Lp{t) = e/h jt^V{t')dt' 
in the phases of the transmission amplitudes, thus making them time dependent. By going 
through the argument presented in Sec. IV, one obtains the scattering states ([T7|) and ([19D 
with time-dependent forward scattering amplitudes: 

ALiR) - AL(ij) e±^^(*^) , (42) 

where the time U = t — \x\/vf is taken retarded in order to account for the finite speed of 
motion after scattering. As before, here we assume that scattering by the potential as well 

on 



as traversing the region where the voltage is apphed takes neghgible time compared to the 
duration of the voltage pulse. In this approximation the amplitudes of backward scattering 
are time- independent constants. 
To draw a relation with the orthogonality catastrophe problem, let us study the effect 
of the voltage pulse on the scattering phases 6i, 62. They can be found by diagonalizing the 
scattering matrix 



(43) 



and writing its eigenvalues as e**^^, e**^^. The relation between the phases 81^2 before and after 
the pulse is written conveniently through 5± = {61 ± 62) /2. The phase 5+ does not change 
at any time, and the phase 6- changes according to 

cos^ 6- {t') + cos^ 5_ (t ) - 2 cos 6- (f) cos 5_ (t) cos Ay^ = | ^ sin^ , (44) 

where Aip = ip{t') — ip{t). Now, let us compare to the orthogonality catastrophe in the Fermi 



system subjected to a time-dependent perturbation (|43|) . Change of the flux induces the 



shift of the phases 6± — > 5^ and makes the new ground state orthogonal to the old one: 

(0'|0)=expf-24lnT^l , (45) 



"7r2 ""A 



F , 



where L is the system size, Xp is Fermi wavelength, and e * is an eigenvalue of the matrix 
S-^{t = oo)S{t = -00): 

sm — = \ Al\ sm — . (46) 

In terms of dynamics, this implies that the old ground state is shaken up so that infinitely 
many particle- hole pairs are excited . It should lead to a logarithmically diverging contri- 



bution to noise, since for each of the particle-hole pairs there is a finite probability (equal to 
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I^L-B/jp) that the particle and the hole will go to different terminals of the conductor, thus 
resulting in a current fluctuation. The periodicity in Faraday's flux $ = — c / V{t)dt follows 
from the gauge invariance and is explicit in Eqs.(0,^) for 6'^. The logarithmic divergence 
vanishes at $ = n^Q, as expected, since at integer $ there is no long-term change of the 
scattering. 

Let us calculate the mean square fluctuation of the charge ((Q^)) transmitted through 



the system due to the pulse. For that, one can use the formula (^61) with the time- dependent 
scattering matrix p3|). To get the second cummulant {{Q^)) one expands the exponent (p6D 
up to second order terms in A, and takes an irreducible average using Wick theorem. The 
averages of Ci^r have the usual form: 



2^ ' 



(47) 



where n{E) = (e^^'^ + 1) ^ is the Fermi distribution. The result reads 



m = 1^/(1^'^ 



+ \AB\ 



^iLp{t')+iujt' 



huj duj 
cucoth— - — 
2T 2-K 



(48) 



where g is spin degeneracy. The first term in (48) is a part of equilibrium noise since it does 
not depend on 93. To analyze the second term, let us take a step-like time dependence of ip 
resulting from an abrupt voltage pulse applied at the time to, < to < ^, the pulse duration 
r being much shorter than t. Taking the integral and keeping only the terms diverging at 
t — > 00, we find 



2-K 



lUO 



lljj 



2, duJ 



\LjJ\ 
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ge^ tEp , ^ . 2 1 ^ 



9 ,ln^ + 2sm' — In- , (49) 

where the ultraviolet- diveregent integrals are cut at frequency ~ EpjTi. By subtracting the 
result for $ = as corresponding to equilibrium, one obtains a logarithmic contribution to 
the non-equilibrium noise: 

((Q^)) = ^AABf \\ sin^ ^ In ^ + + . . . + ((Q))^, , (50) 

The origin of the non-diverging term in Eq. (|50D proportional to $/$o will be discussed 
below. The dots in Eq. (^) represent corrections higher order in $0/"^) and the equilibrium 
noise 

((QX = ^ln^, G = g^|Ar, (51) 

is obtained by repeating the calculation for $ = 0. The expression agrees with the 
Nyquist formula 

(O.J-.)) = e^G uj coth^ (52) 

taken at T = 0, Fourier transformed, and combined with the relation Q = fQj{t')dt'. 

The term in Eq. ( |5D| ) proportional to $/$o is obtained by rewriting the integral in the 
second term of (H) as 

III ^|^|rft^rft2e*^^^*^^-^^*^^+"^*^-*^^\ (53) 

and extracting the contribution of almost coinciding times ti and ^2 by going to new variables 
t = {ti + ^2)72, t' = ti — t2, and changing the order of integrations: 

^1^1 1 ^ 1 1^1^^ ^ (54) 
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where we replaced (p{ti) — (p{t2) = (p{t + t' /2) — (p{t — t' /2) by (p t'. The result ( pi] ) is 
approximate: it does not give the log-term because the transformation (Q) properly takes 
care of the integral only in the domain ti — t2, under the restriction that is varying 
sufficiently smoothly. When is a monotonous function, > 0, the integral in the right 



hand side of (JSJ) equals 27r$/$o and thus produces the term of Eq. (^) proportional to 

It is clear from the derivation that the two terms in the brackets in Eq. (^) arise from 
different integration domains in the ti-t2 space: the first term corresponds to |ti_2| > t, tih < 
0, while the second one is due to almost coinciding moments, \ti — ^2! t. Since these 
domains are almost non-overlapping, the two contributions to the noise (^) do not interfere 
(cross terms are small). 

In order to estimate the correction to the result (|50|), let us derive it by another method 
that allows to trace out the higher order terms. For that, let us take the fiux in the form 
(fit) = NX(t), where X(t) is a smooth monotonous function, A(— 00) = 0, A(oo) = 27r. For 
integer ^ 1 the Fourier component of e*^'^^*-' entering Eq. (|5^) in the stationary phase 
approximation is given by 



-00 ^ ' 



where the dots indicate terms ~ 0(A^~'^/^), and t^'s are real solutions of the equation 

NX{t) + uj = 0. Then we can write 

00 



-00 



and thus obtain 



00 



im . . / E . ... . (57) 



NX{t 



OA 



where the dots represent higher order terms. By differentiating both sides of the equation 
NX{t) = —uj one finds the relation du = —NX{tk)dtk, which means that \uj\duj/X{tk) = 
— \\{tk)\dtk, and therefore the integral in Eq. ( |57D equals N J^^dX = 2txN. Since \uj\dui 
scales as A^^, the correction to Eq. (|57D can be evaluated as 0(1), i.e., it is of the order of 



one for any N. This means that Eq. ([77|) has relative accuracy of 0(1/A^). 

The term in (^) proportional to $/$o is interesting in connection with the picture of 
binomial statistics presented in Sec. VI. In the dc bias case, the distribution of charge for a 
single channel situation was found to be binomial with frequency of attempts equal to eV/ h 
and the probabilities of outcomes p = q = Taken literally, this means that the 

attempts to transfer charge are repeated regularly in time, almost periodic with the period 
h/eV, with each attempt having two outcomes - transmission or reflection - occurring with 
the probabilities p and q. However, the regularity of the attempts does not lead to an ac 
component in the current, rather it appears just as a part of statistical description of charge 
fluctuations. Still, the presence of a non-zero frequency in a non-interacting system requires 
interpretation. 

Let us suppose that the flux varies linearly with time, = —cVt. Since the e.m.f. = 
—d^/cdt, the linear dependence of $(t) is equivalent in its effect on the noise to constant 
voltage V. In accordance with one's expectation, the second term in the brackets in Eq. (|50D 
for a single channel is ((Q^)) = ge^|ylSp$/$0) it is precisely of the form arising from the 
binomial distribution with probabilities of outcomes p and q, and the number of attempts 
= $/$o- (Let us recall that the second moment of the binomial distribution equals pqN.) 
Taking into account that the time during which the flux changes by $o is h/eV, we can 
interpret the number of attempts in the statistical picture as the number of flux quanta by 



which the flux is changed. Such a conclusion suggests an interesting generalization of the 
picture of binomial statistics by attributing the meaning of the number of attempts to the 
flux change measured in the units of $0) regardless of the linear or non-linear character of 
the flux dependence on time. 

It is appealing to put the special role of integer fluxes in connection with the binomial 
statistics of current, where the flux quanta are naturally interpreted as discrete attempts 
to transmit charge. Although this picture is yet to be confirmed by analytic treatment, it 



receives some support from the property of the $0 "Periodic term in ([50| ) to vanish at every 
integer $. One may conjecture that the statistics are close to binomial only when the flux 
change is an integer, and have diverging logarithmic corrections otherwise. The distinction 
that Eq. (^0]) makes between integer and non-integer values of the flux and the relation of 
integer flux change to the number of attempts in the binomial distribution, gives another 
perspective to the statistical picture of a current pulse. 

To summarize, the fluctuations caused by a voltage pulse, in contrast to the average 
transmitted charge, distinguish between integer and non-integer flux change. As a result, 
the dependence of noise on the flux is non-monotonous and has minima at integer values of 
the flux. 

VIII. COHERENT STATES OF CURRENT 

The question we address in this section is about optimal way of changing flux that 
minimizes induced noise. It is clear from what has been said that to achieve minimum of 
the noise one should change the flux by an integer amount. 



9fi 



Ay9 = (p{t = oo) — ip{t = — oo) = 2?™, (58) 

in order to suppress the logarithmically divergent term. However, since for a given A(p the 
noise depends on the actual function f{t), not just on A(f, we have a variational problem to 
solve for the noise as a functional of the time dependence of the flux. This functional was 
derived in Sec. Vll. At zero temperature it is given by 

((g2^) = ^\AB\' f f e^^(*)+-*rft V|— , (59) 
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where A and B are transmission and reflection amplitudes, and g is spin degeneracy. We 
shall study the variational problem (^) with the boundary condition ( |58D and show that 
its general solution has the form of a sum of soliton-like functions: 

m = ±^f: tan-i (^) , r, > 0, (60) 



^ k=i ' ^fc 

where and are arbitrary constants. Under the condition (0), any time dependence of 
the form (|60|) gives absolute minimum to the noise: 



mm 



2\ 



ge^lABl"^ \n\ . (61) 



For an optimal time dependence of the voltage V = —d^/cdt, therefore, one has a sum of 
Lorentzian peaks: 

In order to compare quantum noise with conductance, let us mention that the average 
transmitted charge 

m) = g^l^l'^ = Sll^l'/ ^(t^dt (63) 
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simply obeys the Ohm's law, i.e., there is no particular dependence on the way the flux 
change A<^ is realized. 



of charge fluctuations. For the binomial distribution with probabilities of outcomes p and 
g, p + g = 1, and with the number of attempts N , the second moment is known to be 
equal to pqN . The comparison with Eq. ( |6lD suggests to attribute io n = A$/$o the 
meaning of the number of attempts. This interpretation is supported by the structure of the 
function (^0]) consisting of n terms, each corresponding to unit change of flux. A remarkable 
property of the function (^) is its separability, manifest both in the form of the terms and 
in the way the parameters tfc, enter the expression. Let us note that by making some of 
the tfc's close to each other one can have an overlap in time of different "attempts". The 
overlap, however, does not change the fluctuations (|6TD. The situation reminds the one with 
solitons in integrable non-linear systems, or with non-interacting instantons in integrable 
fleld theories. Also, the absence of interference is interesting in the context of coherent 
nature of transport in this system: after all, we simply have scattering by a time-dependent 
potential. Perhaps, proper interpretation of this effect should be sought in establishing 
relation with the theory of coherent states, known to eliminate to some extent the quantum 
mechanical interference. 

Let us now turn to the variational problem. It is convenient to do the integral over uj 
flrst and to rewrite ( ^9] ) as 



The result (|5DD,(pTD has a simple interpretation in terms of the binomial statistics picture 




(64) 



where D = ^\AB\'^. In order to avoid divergence ai t = t' the denominator in ( |64|) should 
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be understood as 

1 1 



6^0 , (65) 



-{t-t' + i6y {t-t'-i6) 

the condition that one obtains by introducing regularization in (^): \uj\ |co'|e~''^''^. By 
considering variation of the functional (|6^) we have the equation for an extremum: 



Im 



{t - t') 

By using Cauchy formula one checks that the functions 



e^^(*) / -^—l^dt = . (66) 



t-X 



n 7 — ^ ' ^k = tk + irk , (67) 



k=i ^ ~ 



satisfy ( |66D provided that r^'s are all of the same sign. Obviously, the functions (^) are 
just another form of (|60|). 

It remains to be shown that the functional reaches its minimum on the solutions (^^7\j. 
To prove it we proceed in the following steps. Let us write e*'^'-*-' as 

e^^W=/+(t) + /_(t), (68) 

where f+{t) and f~{t) are bounded analytic functions in the upper and lower complex t half- 
plane, respectively. Representation (|68|) exists for any non-singular function and defines the 
functions /+ and /_ up to a constant. Then we substitute Eq. (0) in (§^), and apply 
Cauchy formula for the derivative, 

where the contour of integration is chosen in the half-plane of analyticity of /+ or /_, 
respectively. Thus one gets 

9Q 



-iD iUU -f_f_)dt. (70) 



On the other hand, 

n = — /e-^^wi^^Wcit 
27rz J at 

where the last equahty is a result of substituting (^Sj) and using the relations 

/ /+/- = //-/+ = , (72) 

which follow from Cauchy theorem. Now, Eq. ([7D| ) can be rewritten through Fourier com- 
ponents of /+ and /_ as 

m) = D j {\u{uj)\' + i/-(-^)n ^ ^ , (73) 



thus demonstrating positivity of both terms in (^). (It is used that f+iuj) = /_(— c<j) = 
for u; < 0.) With this, by comparing (^) and ( |7TD we obtain 



((g')) > 27rD \n\ . (74) 

Equality in (^) is reached only when either /+(t) or /-(t) vanishes. Therefore, to obtain 
the minimum one has to take the functions e*'^'-*^ that are regular in one of the half-planes. 
This remark is sufficient to see that the functions ( p^ ) form a complete family of solutions. 

It is worth mentioning that the method used to derive ( |7^ copies almost entirely the 
procedure of derivation of the duality condition in the theory of instantons. Like in other 
situations where the duality condition holds, our "solitons" do not interact: {{Q'^)) shows 
no dependence on the parameters Xk of the solution (|67|) . Among numerous field theories 
that allow for exact solution of the instanton problem the one most similar to our case is 



the theory of classical Heisenberg ferromagnet in two dimensions. For this problem the 
instantons were found by mapping the order parameter space (i.e., the unit sphere) on the 
complex plane The duality condition was shown to take the form of the constraint of 



analyticity or anti-analyticity of the mapped order parameter function (compare with the 
condition /+ = or /_ = derived above). Multi-instanton solutions were given as products 
of single instanton solutions (cf. Eq. (|67|) ). This analogy obviously deserves more attention. 

At this point let us examin an interesting non-optimal time dependence of the flux, the 
sum of two solitons with opposite charge: 

^(t) = 2ftan-V^) -tan-V^)l , (75) 
L \ n ^ r2 

Ti 2 > 0. This function corresponds to e**^*^*^ of the form ( |67D but with the poles in both 
half-planes. In this case Aip = 0, and thus {{Q)) = 0, so min[ ((Q^)) ] = 0. With the function 
([751), however, one finds 

m = 4vrD , (76) 

Ai — A2 

where Ai,2 = ^1,2 + ^^1,2- For different values of the parameters ti^2, ^1,2 Eq. (^) interpolates 
between two trivial limiting cases: 

(i) {{Q"^)) 0, when the two flux steps in ([75|) have nearly the same duration and almost 
overlap; 

(ii) ((Q^)) —>■ AttA, when the fiux steps either differ strongly in their duration or do not 
overlap. 

In the case (ii) the noise is twice bigger than the noise due to a single step, as it should be. 

We see that when Aip/2TT is of the order of one a non-optimal time dependence (p{t) can 
considerably enhance the noise. It is not the case, however, for A(/9/27r ^ 1. This limit was 



studied in Sec. VII, where it was found that when ip{t) is a monotonous function the result 

HQ')) = ge'\AB\'\A^/2n\ (77) 

is rather accurate, even if the time dependence (p(t) is not optimal [^. 

A more intuitive way to understand the accuracy of Eq. ( [77D is to note that for a given 
n the number of parameters in the optimal flux dependence (^) is 2n, which means that 
half of them are in some sense redundant. Because of that any smooth monotonous function 
with sufficiently large variation Aip can be rather accurately approximated by a function of 
the form (|53|), and therefore the noise exceeds the lower bound just slightly. 

An implication of this result for the binomial statistics picture is as follows. As it was 
discussed above there is a (conjectured) correspondence of the terms of Eq. ( |5DD and of the 
attempts. The deviation from the binomial distribution, that of course should exist for a 
non-optimal flux function (p{t), will remain bounded in the case of a smooth (p{t), as Aip 
increases taking integer values. More precisely, the distribution will be written as a mixture 
of binomial distributions with different numbers of attempts, P(m) = J2n PNPN{fn), 
where Pj\f{m) = p"^q^~'^C^. The estimated correction implies that the distribution of 
attempts pn has finite variance in the limit A^ = Aip/2TT — > oo. 

Before closing, let us mention that in order to apply the results of Sees. VII, VIII to 
transport in a mesoscopic metallic conductor with disorder, described by many conducting 
channels with transmission constants T^, one just needs to replace lAB]"^ by I]n^n(l — Tn), 
since different scattering channels contribute to the noise independently. The condition of 
validity of our treatment then is that the variation of the flux is sufficiently slow, so that 
min[rfc] ^ fi/Ec, the time of diffusion across the sample. However, at non-zero tempera- 



ture one also has to satisfy the condition <^ fi/T^ the time of phase breaking. So, the 
temperature interval where our estimate of the noise holds is T < Ec- 

IX. CONCLUSIONS 

We introduced a quantum-mechanical scheme that gives complete statistical description 
of electron transport. It involves a spin 1/2 coupled to the current so that the spin precession 
measures transmitted charge. The off-diagonal part of the spin density matrix, taken as a 
function of the coupling constant, gives the generating function for the electron counting 
statistics. We find the statistics in a single-channel ideal conductor for arbitrary relation be- 
tween temperature and voltage. In equilibrium, the counting statistics are gaussian, both for 
zero-point fluctuations and at finite temperature. At constant voltage and low temperature 
the statistics are BernouUian and the distribution is binomial. 

The theory leads to interesting conclusions applied to the current fluctuations produced 
by a voltage pulse. In this case, the noise has phase sensitivity: it oscillates as function 
of Faraday's flux, c jV{t)dt, reaching minimum at integer fluxes. We studied the noise 
as function of the shape of the voltage pulse and found optimal time dependence that 
provides absolute minimum of the noise for given average transmitted charge. Solution 
displays interesting analogy with the problem of instantons in the field theories obeying 
duality symmetry. Optimal time dependence is a sum of Lorentzian peaks of voltage, each 
corresponing to a sohton of flux. The change of flux for a soliton is equal to the flux 
quantum $o- The solitons are interpreted in terms of the binomial statistics picture of 
charge fluctuations as attempts to transmit electrons, one electron per soliton. 



APPENDIX A: LARMOR CLOCK MEASUREMENT OF TUNNELING TIME 



How long does it take a particle to tunnel under a barrier? More precisely, suppose a 
particle of energy E is moving in one dimension, and is scatterred on a potential barrier: 
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2 dx' 



+ U{x) 



ilj{x,t). (Al) 



What is the probability that during the scattering the particle spends time r within the 
region a < x < b under the barrier? Questions of that kind arise naturally in discussion of 
any quantum-mechanical process that takes finite time, like nuclear or chemical reactions, 
resonance scattering, or tunneling. 



There have been several attempts to treat such problems that resulted in formula- 
tion of a very interesting concept of Larmor clock. It has various analogies with the spin 
galvanometer discussed above, and it seems useful to review the Larmor clock here using 
the same language. The Larmor clock uses an auxiliary spin 1/2 attached to the scattering 
particle, and an auxiliary constant magnetic field a; localized within the region of interest, 
a < X < b, 

1 } 

'Hint = ~2^^^j i^^{x)i){x)dx . (A2) 

a 

The choice of coupling is such that the spin precession angle is proportional to the time 
spent in the region a < x < b. The difference from our spin-galvanometer is that the spin is 
not stationary, but travels with the particle, and also that the spin is coupled to the particle 
density, rather than to the current. 

To find the distribution of times one has to write down the system density matrix evolved 
in time, and take partial trace over the particle outgoing states. (We assume that one 



does not have to distinguish between different results of scattering, and is interested in the 
tunnehng time only, regardless of whether the particle went through the barrier, or has been 
reflected.) Then, by following the argument of Sec. Ill one obtains the spin density matrix: 



Psit) 



PTt(O) x(^)PTi(0) 



(A3) 



Here 



(A4) 



where e is the evolution operator for the one-particle problem with no spin: 



ij{x,t) , 



(A5) 



where 9ab = 0{x — a)6{h — x). The auxiliary magnetic field ui now turns into a constant 
potential within the region a < x < h. Here again, with the spin degrees of freedom taken 



care of by 



trace one finds 



we are left with a single particle problem. By using cyclic property of the 



X(^) = (e 



(A6) 



Here the brackets (...) mean averaging over the particle initial state. Note that xil^) is 
written in terms of a purely single particle problem, not involving spin variables. 

The quantity xi}^) obtained by measuring precession of the spin is a generating function 



for the distribution of times, which is clear from the Fourier transform 



x(^) 



(A7) 



The probabilities -P(r) of different precession angles of the spin should be interpreted as the 
scattering time distribution. 



The probabilities -P(t) defined by ( |A5| ) , 
ties: 

a) JP{T)dT = l; 

b) P{t) are real numbers; 

c) P{t) vanish at negative times r < 0. 
The normalization property is derived from by setting uj = 0. Property b) (real 
valuedness) is derived from x(— cj) = x{^) which follows from (|A6|) . The causality property 
c) follows from considering the evolution in the problem with u continued to complex 
values. One notes that both the solution ip{x,t) of Eq. (^) and the evolution operator 
^-iHuit g^j,g regular in the upper half-plane Im > 0, which means that the same is true 
for x(co'). From that, the causality property c) follows by the usual argument using Cauchy 
theorem in the integral 

P(r) = X(c^)e-^-^^ (A8) 

by closing the integration contour in the upper half-plane. 

The properties a), b) and c) suggest that P{t), so far defined formally as Fourier spec- 
trum of x{^)j can have a meaning of probability. However, generally the sign of P{t) can 
be either positive or negative, which makes the probabilistic interpretation problematic. 

For the one particle problem one can write the generating function x{^) terms of the 
scattering amplitudes A and B. For that, it is convenient to use the expressions (p3D, 
for the evolution operator in terms of the scattering matrix S, written using the wave-packet 
scattering states (|T^. Specializing to one particle and taking partial trace, one finds 

Xioo) = A^^{E)A^{E) + B_^{E)B^{E) , (A9) 



), and (^1) have several interesting proper- 



where A{u) and B{uj) are the transmission and reflection amphtudes of the problem ( |A5|) 
taken at the energy E of incident particle. 

To see the Larmor clock working, let us consider an example of resonance scattering, 
where a particle is scattered on a potential forming a quasibound state of life-time F. Using 
the method described above one can find the distribution of times it takes the particle to 
scatter. For simplicity, suppose that the particle can be only reflected, but not transmitted 
{A = 0). Then the reflection amplitude as function of energy is given by the Breit-Wigner 
formula: 

Turning on the field u in the quasibound state region is equivalent to shifting the resonance 
energy: Eq Eq — uj/2. Thus, the generating function of the time distribution is 

e — uj — iL e + uj + tl 
where e = 2{E — Eq). The distribution P(t) is found by Fourier transform: 

.duj 
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e 

d . 4F 



P(r) = / x(^)e-'-V 



5{t) (F sin er — e cos er) e 



-Ft 

^ sinere"'^^) . (A12) 

or \ e J 

The 5-function term corresponds to the non-resonance scattering channel. Other terms 
describe dwelling in the quasibound state. In this example P{t) is changing sign, which 
makes the probabilistic interpretation ambiguous. 

The paradox arising due to negative -P(t) is only an apparent one. Really, the measure- 
ment of time performed by the Larmor clock is not the usual quantum-mechanical measure- 
ment, since the time is not an operator, and thus it cannot be measured in the same sense 
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as other quantum-mechanical observables. This should be contrasted with the measurement 
of charge described above. Although the spin precession measurement scheme we use looks 
quite similar to the Larmor clock, there is a difference: Electric charge is an observable in 
the usual quantum-mechanical sense, it takes quantized integer values, and the probabilities 
of those values resulting from our calculation are non-negative. 

APPENDIX B: BOSONIZATION CALCULATION OF COUNTING STATISTICS 

In order to find generating function of counting statistics for a single channel conductor, 
we have to evaluate 

x{\) = {expi\N,) , (Bl) 

where Nt = /(ci^^Ci^^C2^^C2,r)'^T, and q^t-, cf^^ are canonical Fermi operators. 

In one dimension, there is an equivalence between ideal Fermi gas and harmonic Bose 
chain, which provides a representation of the Fermi problem in terms of free bosons, known as 



the bosonization transformation [24-pq]. This representation facilitates calculting averages 



like (pT|), since they are being transformed to the form of a gaussian average . 

According to the bosonization theory, bosonic Hamiltonian representing the fermionic 
problem is written as 

HBose = ^J-- Ml? ■■ + : i^^e^f : dx , (B2) 
where 0^fj^{x) are Bose operators, 

[VeLiR^{x),eLiii){y)] = ±2m6{x-y) . (B3) 



Connection to the fermionic problem is given as a relation between the densities of the left- 



and right-moving fermions, pi{x) = cf^Ci^x, = 1,2, and the bosonic variables OL(^Ry 



X) 



written as 



(B4) 

One notes that the operator A^^ in (|El| ) is linear in the densities pi, and thus it is represented 
by an expression linear in the bosonic variables, 



ZTT ^ 



(B5) 



which turns the average in ( |B1[ ) into a gaussian type. 

Therefore, the average of (^8]) is equal to the product of averages 

X(A) = (exp^(^z.(t) -^L(0)))(exp^(^H(t) -^h(O))) 



(B6) 



taken over the ground state of the Hamiltonian (P2|). To perform the average in ([B6|), it 
is sufficient to deal with the average over ^l's, because of the left-right symmetry of the 
problem. 

Let us write 9l{x) in terms of bosonic operators of plane waves: 



fc>0 



k>0 



(B7) 



One checks that the commutation relations ( p3|) are consistent with canonical commutation 
relations between bk and 6^/. (The Hamiltonian of left-moving fermions is represented by 
T^L = Y.k>o vkb'^bk.) The quantity OL{t) — Ol{0) appearing in the average (p^ ) is written as 



fc>0 



1% + ie-''"" - l)bt 



We evaluate the average 



(((^lW - ^l(0))2)) = E ^le*'^^* - l\\2N^Ukv/T) + 1) 



fe>0 



/UK /I \ 

— sin2(wH/2) coth(i;A;/2T) = 2 In — - sinh(7rTt) . (B9) 
\k\ \nTo J 

—ao 

This expression equals (27r)^ times the function f{t,T) computed in (^). From that, we 
find the average ( P^ ) to be 

Ztt 

= exp[-{^ymit)-e,{0)Y))] =exp(-AV(t,T)) , (BIO) 

which is the desired result. 

Periodicity of xW corresponding to the charge quantization, is recovered if one 

corrects the relation between pi{x) and 6i{x), in order to take into account the integer- 
valuedness of the particle number A^^. Using the relation, pi{x) = and performing 

n 

the average, one arrives at the result (p9|). 
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